An Effective Membrane Model of the Immunological Synapse 
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The immunological synapse is a patterned collection of diflFerent types of receptors and ligands 
that forms in the intercellular junction between T Cells and antigen presenting cells (APCs) during 
recognition. The synapse is implicated in information transfer between cells, and is characterized 
by different spatial patterns of receptors at different stages in the life cycle of T cells. We obtain 
a minimalist model that captures this experimentally observed phenomenology. A functional RG 
analysis provides further insights. 

PACS numbers: 87.16.Dg,64.60.-i 



T lymphocytes (T cells) are the orchestrators of the adaptive immune response in complex organisms. A key event 
during activation of the immune response is T cell recognition of cells that display peptides derived from foreign 
antigens on their surface . Recent experiments H, ^ IS have vividly demonstrated that during this recognition 
process a highly organized pattern of different types of receptors and ligands forms in the intercellular junction between 
T cells and antigen presenting cells. This recognition motif is several microns in diameter, and since it is implicated 
in information transfer between the cells, it is called the immunological synapse. Formation of a synapse is also 
characteristic of an earlier stage in the life cycle of T cells. Immature T cells (thymocytes) are selected in the thymus 
so that they are not activated by peptides derived from the organism itself 0|- In the thymus, thymocytes interact 
with cells that display self peptides on their surface. Thymocytes that bind strongly are deleted by apoptosis. The 
synapses formed during thymocyte selection are distinctly different in character U Q from those observed during 
mature T cell activation. Understanding the mechanisms via which synapses form under different circumstances and 
the biological purpose of creating different spatial patterns of cell surface receptors are active areas of research. In 
addition to the biological significance, an understanding of these issues may also inspire the creation of synthetic 
mimics that could carry out biomimetic recognition tasks which could be useful in applications such as targeted drug 
delivery. 

In this letter, starting from a model proposed by Qi et al. Q, we develop a minimalist model that captures some of 
the essential physics of synapse formation when apposing membranes contain complementary pairs of receptors and 
ligands. The model allows us to calculate a phase diagram which delineates the conditions that lead to a transition 
from synaptic patterns characteristic of mature T cells to those observed during thymocyte selection. This phase 
diagram may serve as a guide for the design of synthetic analogs. 

Consider two membranes containing complementary pairs of receptors and ligands. The intramembrane motion of 
receptors and ligands is determined either by diffusion or a directed velocity toward the center of the junction [lol |. 
Complementary receptors and ligands can bind to each other if apposed. Different receptor-ligand complexes have 
different topographical size [Tlj . and hence the rate of association depends upon local intermembrane separation. 
The receptor-ligand complexes also dissociate at a prescribed rate. For the case where there are only two types of 
receptors- ligand pairs (TCR-MHCp and LFAl-ICAMl in the case of the T cell immunological synapse), the pertinent 
equations are: 
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Here, C represents the concentration of a given molecule, and z is the separation of the T cell membrane from the 
APC membrane. The abbreviations used for different protein molecules are: TCR T, MHCp M, TCR-MHCp 
TM, LFAl ^ L, ICAMl -> I, and LFAl-ICAMl LI. The different species diffuse with the corresponding diffusion 
coefficients D. V is a directed velocity of TCR due to cytoskeletal motion ^Td\. Since small values of V do not change 
the qualitative physics [T^ . for simplicity we do not consider this further in this letter. The binding and dissociation 
chemical reactions are controlled by the rate constants kon{z) and k^ff. The local intermembrane separation, z, 
evolves according to a time-dependent Ginzburg-Landau equation, 

where M is a phenomenological parameter for the rate at which membrane shapes respond to changes in the free 
energy, F = J d?x [^(Vz)^ + ^(V^z)^ + CtmBtm{z) + CliBli{z)] . The first two terms are due to costs associated 
with deforming the membranes: If one of the membranes is supported and undeformed, 7 and k are the surface tension 
and bending rigidity of the free membrane; for similar membranes they are half the corresponding parameters for each 
membrane. Btm{z) and Bli{z) denote respectively the bond energy gained from forming TM and LI complexes for 
a specific intermembrane separation z. 

The time scales associated with large-scale membrane shape changes can be much slower than those associated with 
binding kinetics and protein mobility. In such cases, for a given membrane shape, the local species concentrations are 
related by 
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Substituting Eq. Q in Eq. (j^J obtains the following effective dynamical equation for z 
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This corresponds to model A dynamics (Eq. ^) for the order parameter z with an effective free energy functional 
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The reduction of the full set of equations to one determining the separation z is reminiscent of the elimination of 
counter-ions in constructing the Poisson-Boltzmann equation, where the concentration of counter-ions is obtained 
from the Boltzmann weight of a potential, itself determined self-consistently from the counterion charge profile. Here 
z plays the role of the potential, while the protein complexes are analogous to the counterions. 

Bound receptor-ligand complexes do not favor membrane shape changes that deform the bond from its natural 
length. We thus use a harmonic approximation for the functional form of the bond elasticity energy, i.e. Bi{z) = 
B^ -f \i{z — Zi)^/2, with i = 1, 2 for the two possible complexes 9]. The maximal bond energies B^ can be absorbed 
into Kon /Koff, and shall be ignored henceforth. The natural lengths of receptor-ligand proteins are taken to be those 
appropriate for the T cells (zi = ztm = 15nm and Z2 = zli = 42nm) [I^jIIjI, with corresponding 'spring constants' 
of Xi given below. With this approximation, the dynamical Eq. (4) takes the form (measuring time in units such that 
M = 1) 
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with an effective potential 

^(,).,,^|-^i*exp{-l£^}-C|exp{-(£^}, for.>0 
I c, for z < 0. 

In writing Eq. |(7|) we have assumed that — ksT /\, which is appropriate if the curvature of the receptor-hgand 
interaction potential is the same in the bound state and at the activation barrier. Other assumptions do not alter the 
qualitative physics. The parameters CI = {Kon/ Kof f)CTCM and C| = {K'^^/ K'^^ j)ClCi represent the number of 
complexes (TM and LI, respectively) formed if the intermembrane separation equals the corresponding optimal bond 
separation. As depicted in Fig. 1, the effective potential has two minima corresponding to the two possible bonds. 
Negative separations are prevented by a large energy cost fc^Tc ^ 1 (hard wall in simulations). 

Eq. © was solved numerically by a finite difference method. Initially, the upper membrane was assigned a parabolic 
shape in the contact region, and the intermembrane distance was kept constant (^ 50nm) outside the contact region. 
For the membrane parameters we use experimentally determined values |14| of 7 = lOOksT / ^m? and k — AQQksT^ 
and for the bond elasticity we take cti = limn and (T2 — 5nm 9\. Simulations were carried out for various ratios 
of C*/C|, and local concentrations of the TCR-MHCp and LFAl-ICAMl were obtained from the intermembrane 
separation using Eqs. 

Zero Noise: In the absence of noise ((^ = 0), Eq. © is a relaxational dynamics toward the minimum of the potential 
V{z), subject to the imposed boundary conditions. For > C2 the minimum at ztm is deeper, and indeed in the 
numerical simulations a central cluster of TCR-MHCp complex forms quickly in the intercellular contact area. As 
the boundary of the contact region is kept at a constant height close to zli, the longer LFAl-ICAMl moves toward 
the boundary giving rise to a structure very similar to the final synaptic pattern obtained in the experimental and 
numerical studies 0, 0, 0, ■ 

When the concentration C* is below C|, the global minimum of V{z) switches to z^j. In this regime, the receptor- 
ligand complex LFAl-ICAMl concentration dominates the entire intercellular junction, and no synaptic pattern is 
formed. In the noise-less (mean-field) limit the location of the transition between the two patterns is always at 
C1/C2 = 1. The patterns obtained and the conditions for transitions between them are quite different when C is 
finite. 

Finite Noise: In this case, numerical simulations of Eq. ^ are performed with a random C, with zero mean, and 
variance equal to fc^T to mimic thermal noise. Starting with sufficiently large values of Cj* > C|, we again find 
that TCR-MHCp aggregates at the center of the intercellular junction at long times. In Fig. (2a) a cross section of 
the upper membrane and the associated protein concentrations are shown. Very similar synaptic patterns have been 
observed for in vivo and in vitro experiments with mature T cells. 

For double positive (DP) thymocytes, the TCR concentration is much lower compared to mature T-cells This 
experimental fact corresponds to lower values of CJ. When we numerically solve Eq. © with CI/C2 < fc, the mature 
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FIG. 1: The effective potential V{z) in Eq. Q. 
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FIG. 2: Right panel: (a) Structure of a mature synapse: TCR-MHCp occupies the central region (represented in green), 
surrounded by a ring of LFAl-ICAMl (shown in red), (b) Dynamic accumulation of TCR-MHCp (green) at random places in 
a sea of LFAl-ICAMl (red) molecules as observed for thymocytes. Left panel: Height fluctuations corresponding to (a) and 
(b). The green {ztm = 15nm) and red {zli = 42nm) lines are meant to guide the eye. 

synapse is never formed. Rather, we observe (see Fig. 2b) fluctuating patterns with fleeting clusters of TCR-MHCp 
forming at various locations in a sea of LFAl-ICAMl. Note that this pattern is different in character from the one 
observed for C'l /Cn < 1 in the noiseless case. Such dynamic patterns have been observed in recent experiments using 
DP thymocytes [JQ, and can be understood using the simple picture of relaxational motion in a potential with 
two minima. As we lower the TCR concentration, the minimum at z — zxm becomes shallower, and eventually the 
minimum near z — zn corresponds to the global equilibrium. However, thermal fluctuations allow sampling of the 
region near the minimum a.i z — ztm- This will favor sporadic aggregation of TCR-MHCp at random places. As long 
as the barrier between the two minima is not very large, such transient fluctuations will be frequent. 

The effects of thermal fluctuations can be estimated by expanding V{z) to quadratic order around each of its 
minima. The total energy cost of a deformation of wavenumberp is given by Ei{p) = {kBTC*(7~^ -f -I- Kp'*) |z(p)p. 
Integrating out these (Gaussian) fluctuations yields a free energy per unit area of 



The location of the transition between the two classes of patterns can now be estimated from F* = F2 ■ Thermal noise 
enhances the minimum with larger a (width), because of its larger entropy. Hence the more flexible bond is formed 
more easily when thermal noise is present |l7| . Fig. 3 shows that the larger the temperature, the smaller the value of 
fc- This "phase diagram" could serve as a guide for the design of synthetic systems as it shows how to manipulate 
conditions such that synaptic patterns characteristic of mature T cells or thymocytes are realized. 

Self consistency of the above quadratic expansion requires that the extent of the fluctuations should be less than 
the corresponding variance of the bond. The mean square width of the interface trapped around z = Zi \s given by 



Typical values of CI range from O(10)(/im)^^ for thymocytes to O(100)(/im)^^ in mature T cells. For this range of 
parameters, the constant term in the denominator of Eq. dominates, until it is cut-off by the bending rigidity at 
distances of the order of l/im. (The overall size of the synapse is of the order of lO/im Q.) It is then easy to estimate 
the value of the integral, and we find that the ratios Wi/ai and W2/cr2 are of the order of 0.1. While these estimated 
fiuctuations are smaller than the width of the attractive potentials, they are significant enough to make nonlinear 
corrections to the above quadratic expansion important. Thus the phase diagram in Fig. 3, using Eq. (jS)) should be 
regarded as an approximation. 
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FIG. 3: Phase diagram (projection on the Ci-CJ plane) for physiological temperature T — 37°C and room temperature 
T — 27°C. Ci and CI are measured in units of {fim)~'^. The dotted portion of the curves (for low values of Ci-CJ) is 
conjectured as explained in the text. 



When anharmonic corrections are significant, an alternative approach is to employ a functional renormalization 
group (RG) scheme used 0,0,1131 in the context of wetting transitions. It consists of the following steps: (1) 
The fluctuating field z{x) is divided into two parts corresponding to small wavcnumbers (0 < \p\ < A/b) and large 
wavenumbers (A/6 < \p\ < A); (2) The large wavenumber fluctuations A/b < |p| < A are integrated out to yield a 
coarse grained Hamiltonian; and (3) The system is rescaled to new coordinates x' — x/b and z'{x') = z{x). (The field 
z does not need to be scaled in two dimensions.) The second step can only be done approximately, and a first order 
expansion of the potential V{z) leads to the linear approximation^^ 
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Note that k is an irrelevant operator in the RG sense (vanishing as b~^ under scaling). Indeed the second term in the 
above equation can be ignored asymptotically for 6^1, but may be important for quantitative comparisons when 
the bending rigidity is large. 

Applying Eq. (|10|l to the potential in Eq. {Tj) yields the renormalized potential 



Vb{z) « b^ksT 

The soft wall renormalizes into an error function of width S, while the variance of the Gaussian potential is increased 
by 6^. Thus, the bigger the width of the potential, the deeper it gets due to the fiuctuations, as it allows more 
flexibility (entropy) in bond formation. For large values of the width S (highly fluctuating surfaces), the location of 
the transition point is shifted to CJ'/CI = fc = o'2/(Ti ~ 0.38. Due to the factor of b^ in Eq. (|12|l the net potential 
gets stronger under RG, and the linear approximation eventually breaks down. 

The RG transformation should be extended to a rescaling factor b at which the curvature at the minimum of the 
renormalized potential is such that the resulting fluctuations in width are comparatively small. Beyond this scale, 
a quadratic expansion around the minimum of the effective potential should be valid. As demonstrated previously, 
typical parameters for the synapse yield a bare potential at which the above criterion is (barely) satisfled. Thus the 
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linear RG approximation is of limited practical value in this context. It may, however, yield interesting information 
about the global form of the phase diagram as the concentrations Cj and are further reduced. From the double 
well shape of the potential, we anticipate that the overall phase diagram consists of a line of discontinuous transitions 
separating preferences in the synapse for TM or LI bonds. This line is located at Cj^/Ca < 1 for large concentrations, 
and moves to C^/CI « axja-i at lower concentrations. It presumably terminates at an Ising critical point when the 
fluctuations in width become of the order of the separation ztm — zli between the two minima. Sufficiently small 
values of CJ' and C| should also result in an unbinding transition which is located at Clai + C|c72 oc \fk^Tl^. 
It is thus useful to make a more detailed study of the phase diagram of the model synapse as well as the dynamic 
trajectories as various phase boundaries are approached. 
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